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Abstract 

A grove is a spanning forest of a planar graph in which every component tree con- 
tains at least one of a special subset of vertices on the outer face called nodes. For the 
natural probability measure on groves, we compute various connection probabilities for 
the nodes in a random grove. In particular, for "tripartite" pairings of the nodes, the 
probability can be computed as a Pfafnan in the entries of the Dirichlet-to-Neumann 
matrix (discrete Hilbert transform) of the graph. These formulas generalize the deter- 
minant formulas given by Curtis, Ingerman, and Morrow, and by Fomin, for parallel 
pairings. These Pfafnan formulas are used to give exact expressions for reconstruction: 
reconstructing the conductances of a planar graph from boundary measurements. We 
prove similar theorems for the double-dimer model on bipartite planar graphs. 



1 Introduction 

In a companion paper [KW06J we studied two probability models on finite planar graphs: 
groves and the double-dimer model. 



1.1 Groves 

Given a finite planar graph and a set of vertices on the outer face, referred to as nodes, a 
grove is a spanning forest in which every component tree contains at least one of the nodes. 
A grove defines a partition of the nodes: two nodes are in the same part if and only if they 



are in the same component tree of the grove. See |Figure 1 

When the edges of the graph are weighted, one defines a probability measure on groves, 
where the probability of a grove is proportional to the product of its edge weights. We proved 
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Figure 1: A random grove (left) of a rectangular grid with 8 nodes on the outer face. In 
this grove there are 4 trees (each colored differently), and the partition of the nodes is 
{{1}, {2, 7, 8}, {3, 4, 5}, {6}}, which we write as 1 1 278 1 345 1 6, and illustrate schematically as 
shown on the right. 



in |KW06| that the connection probabilities — the partition of nodes determined by a random 
grove — could be computed in terms of certain "boundary" measurements. Explicitly, one 
can think of the graph as a resistor network in which the edge weights are conductances. 
Suppose the nodes are numbered in counterclockwise order. The L matrix, or Dirichlet-to- 
Neumann matrij0 (also known as the response matrix or discrete Hilbert transform), is then 
the function L = (L^j) indexed by the nodes, with Lv being the vector of net currents out of 
the nodes when v is a vector of potentials applied to the nodes (and no current loss occurs 
at the internal vertices). For any partition tt of the nodes, the probability that a random 
grove has partition n is 

Pr(7r) =Pr(7r)Pr(l|2|---|n), 

where 1|2| • • • |n is the partition which connects no nodes, and Pr(7r) is a polynomial in 
the entries Li j with integer coefficients (we think of it as a normalized probability, Pr(7r) = 
Pr(7r)/ Pr(l|2| • • • \n), hence the notation). In [KW06] we showed how the polynomials Pr(7r) 
could be constructed explicitly as integer linear combinations of elementary polynomials. 

For certain partitions ir, however, there is a simpler formula for Pr(-7r): for example, 
Curtis, Ingerman, and Morrow [CIM98] , and Fomin [FomOlj . showed that for certain parti- 
tions 7T, Pr(7r) is a determinant of a submatrix of L. We generalize these results in several 
ways. 



Firstly, we give an interpretation ( § 8 ) of every minor of L in terms of grove probabilities. 



This is analogous to the all-minors matrix-tree theorem [Cha82j |Che76[ pg. 313 Ex. 4.12 



4.16, pg. 295], except that the matrix entries are entries of the response matrix rather than 
edge weights, so in fact the all-minors matrix-tree theorem is a special case. 



1 Our L matrix is the negative of the Dirichlet-to-Neumann matrix of |CdV98j . 
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Secondly, we consider the case of tripartite partitions it (see Figure 2 ), showing that 
the grove probabilities Pr(7r) can be written as the Pfaffian of an antisymmetric matrix 
derived from the L matrix. One motivation for studying tripartite partitions is the work 
of Carroll and Speyer ( 'SO 1 and Petersen and Speyer [PS05] on so-called Carroll- Speyer 
groves (Figure 7) which arose in their study of the cube recurrence. Our tripartite groves 
directly generalize theirs. See 



9 



A third motivation is the conductance reconstruction problem. Under what circum- 
stances does the response matrix (L matrix), which is a function of boundary measure- 
ments, determine the conductances on the underlying graph? This question was studied in 
[UIM98L ICdV98l ICdVCV96j . Necessary and sufficient conditions are given in [CdVCV96j 
for two planar graphs on n nodes to have the same response matrix. In |CdV98] it was 
shown which matrices arise as response matrices of planar graphs. Given a response ma- 
trix L satisfying the necessary conditions, in § 7 we use the tripartite grove probabilities to 



give explicit formulas for the conductances on a standard graph whose response matrix is L. 
This question was first solved in |CIM98] , who gave an algorithm for recursively computing 
the conductances, and was studied further in |CM02t IRusQ3j . In contrast, our formulas are 
explicit. 




Figure 2: Illustration of tripartite partitions. The two partitions in each column are duals 
of one another. The nodes come in three colors, red, green, and blue, which are arranged 
contiguously on the outer face; a node may be split between two colors if it occurs at the 
transition between these colors. Assuming the number of nodes of each color (where split 
nodes count as half) satisfies the triangle inequality, there is a unique noncrossing partition 
with a minimal number of parts in which no part contains nodes of the same color. This 
partition is called the tripartite partition, and is essentially a pairing, except that there 
may be singleton nodes (where the colors transition), and there may be a (unique) part of 
size three. If there is a part of size three, we call the partition a tripod. If one of the 
color classes is empty (or the triangle inequality is tight), then the partition is the "parallel 
crossing" studied in [CIM98] and |Fom01] . 
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1.2 Double-dimer model 

A number of these results extend to another probability model, the double-dimer model on 
bipartite planar graphs, also discussed in [KW06J. 

Let Q be a finite bipartite graph! embedded in the plane with a set N of 2n distinguished 
vertices (referred to as nodes) which are on the outer face of Q and numbered in counter- 
clockwise order. One can consider a multiset (a subset with multiplicities) of the edges of Q 
with the property that each vertex except the nodes is the endpoint of exactly two edges, and 
the nodes are endpoints of exactly one edge in the multiset. In other words, it is a subgraph 
of degree 2 at the internal vertices, degree 1 at the nodes, except for possibly having some 
doubled edges. Such a configuration is called a double-dimer configuration; it will connect 
the nodes in pairs. 

If edges of Q are weighted with positive real weights, one defines a probability measure 
in which the probability of a configuration is a constant times the product of weights of 
its edges (and doubled edges are counted twice), times 2 l where I is the number of loops 
(doubled edges do not count as loops). 

We proved in [KW06j that the connection probabilities — the matching of nodes deter- 
mined by a random configuration — could be computed in terms of certain boundary mea- 
surements. 

Let Z DD ((?,N) be the weighted sum of all double-dimer configurations. Let Q BSN be the 
subgraph of Q formed by deleting the nodes except the ones that are black and odd or white 
and even, and let Qff be defined as Q BW was, but with nodes i and j included in Qff if 
and only if they were not included in Q . 

A dimer cover, or perfect matching, of a graph is a set of edges whose endpoints cover 
the vertices exactly once. When the graph has weighted edges, the weight of a dimer con- 
figuration is the product of its edge weights. Let Z BSN and Zfj be the weighted sum of 
dimer configurations of Q BVJ an Qff , respectively, and define Z WB and Zf B similarly but 
with the roles of black and white reversed. Each of these quantities can be computed via 
determinants, see |Kas67j . 

One can easily show that Z DD = z BW Z WB ; this is essentially equivalent to Ciucu's graph 
factorization theorem [Ciu97j . (The two dimer configurations in Figure 3 are on the graphs 



G BW and Q SNB .) The variables that play the role of Ljj in groves are defined by 

Xij = Z B f/Z BW . 



We showed in |KW06j that for each matching tt, the normalized probability Pr(7r) = 
Pr(7r)Z WB /Z BW that a random double-dimer configuration connects nodes in the matching tt, 
is an integer polynomial in the quantities X^j. 

2 Bipartite means that the vertices can be colored black and white such that adjacent vertices have different 
colors. 



Combinatorics of tripartite boundary connections for trees and dimers 



Kenyon & Wilson 



5 




In the present paper, we show in ITheorem 6.11 that when ix is a tripartite pairing, that 
is, the nodes are divided into three consecutive intervals around the boundary and no node 
is paired with a node in the same interval, Pr(7r) is a determinant of a matrix whose entries 
are the Xi/s or 0. 



1.3 Conductance reconstruction 

Recall that an electrical transformation of a resistor network is a local rearrangement of the 



type shown in Figure 4 These transformations do not change the response matrix of the 
graph. [CdVGV96] showed that a planar graph with n nodes can be reduced, using electrical 
transformations, to a standard graph S n (shown in |Figure 5| for n up to 6), or a minor of 
one of these graphs (obtained from S n by deletion/contraction of edges). 

In particular the response matrix of any planar graph on n nodes is the same as that for 
a minor of the standard graph E n (with certain conductances). |CdV98] computed which 
matrices occur as response matrices of a planar graph. [CIM98 ] showed how to reconstruct 




a 



a 
+ 





a+b+c 



Figure 4: Local graph transformations that preserve the electrical response matrix of the 
graph; the edge weights are the conductances. These transformations also preserve the 
connection probabilities of random groves, though some of the transformations scale the 
weighted sum of groves. Any planar graph with n nodes can be transformed to a minor of 
the "standard graph" S n ( Figure 5 ) via these transformations [CdVG V96] . 
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Figure 5: Standard graphs S n with n nodes for n up to 6. 



recursively the edge conductances of S n from the response matrix, and the reconstruction 
problem was also studied in [CM0 2] and |Rus03j . Here we give an explicit formula for the 
conductances as ratios of Pfaffians of matrices derived from the L matrix and its inverse. 
These Pfaffians are irreducible polynomials in the matrix entries (ITheorem 5.11) . so this is in 
some sense the minimal expression for the conductances in terms of the Lij. 



2 Background 

Here we collect the relevant facts from |KW06j . 

2.1 Partitions 

We assume that the nodes are labeled 1 through n counterclockwise around the boundary 
of the graph Q. We denote a partition of the nodes by the sequences of connected nodes, 
for example 1 1 234 denotes the partition consisting of the parts {1} and {2,3,4}, i.e., where 
nodes 2, 3, and 4 are connected to each other but not to node 1. A partition is crossing if 
it contains four items a < b < c < d such that a and c are in the same part, b and d are 
in the same part, and these two parts are different. A partition is planar if and only if it is 
non-crossing, that is, it can be represented by arranging the items in order on a circle, and 
placing a disjoint collection of connected sets in the disk such that items are in the same 
part of the partition when they are in the same connected set. For example 13 1 24 is the only 
non-planar partition on 4 nodes. 

2.2 Bilinear form and projection matrix 

Let W n be the vector space consisting of formal linear combinations of partitions of {1, 2, . . . , n}. 
Let U n C W n be the subspace consisting of formal linear combinations of planar partitions. 

On W n we define a bilinear form: if r and a are partitions, (r, cr) t takes value 1 or and 
is equal to 1 if and only if the following two conditions are satisfied: 

1. The number of parts of r and a add up to n + 1. 

2. The transitive closure of the relation on the nodes defined by the union of r and a has 
a single equivalence class. 
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For example (123|4, 24|l|3) t = 1 but (12|34, 12|3|4) t = 0. (We write the subscript t to 
distinguish this form from ones that arise in the double-dimer model in §6 ) 

This form, restricted to the subspace U n , is essentially the "meander matrix", see [KW06[ 
IDFGG97] . and has non-zero determinant. Hence the bilinear form is non-degenerate on U n . 
We showed in |KW06] . Proposition 2.6, that W n is the direct sum of U n and a subspace K n 
on which (,) t is identically zero. In other words, the rank of (,) t is the n th Catalan num- 
ber C n , which is the dimension of U n . Projection to U n along the kernel K n associates to 
each partition r a linear combination of planar partitions. The matrix of this projection is 
called V®. It has integer entries [KW06J. Observe that V® preserves the number of parts 
of a partition: each non-planar partition with k parts projects to a linear combination of 
planar partitions with k parts (this follows from condition 1 above). 

2.3 Equivalences 

The rows of the projection matrix V® determine the crossing probabilities, see lTheorem 2.51 
below. In this section we give tools for computing columns of P^. 

We say two elements of W n are equivalent (=) if their difference is in K n , that is, their 
inner product with any partition is equal. We have, for example, 

Lemma 2.1. 1(234 + 2|134 + 3|124 + 4|123 = 12|34 + 13(24 + 14|23 

which is another way of saying that 

P (t) (13|24) = 1(234 + 2(134 + 3(124 + 4(123 - 12(34 - 14(23. 

This lemma, together with the following two equivalences, will allow us to write any 
partition as an equivalent sum of planar partitions. That is, it allows us to compute the 
columns of V®. 

Lemma 2.2. Suppose n > 2, r is a partition of 1, . . . , n — 1, and r = J2 a ■ Then 

r\n = a a a\n. 

a 

If r is a partition of 1, . . . , n — 1, we can insert n into the part containing item j to get 
a partition of 1, . . . , n. 

Lemma 2.3. Suppose n > 2, r is a partition of 1, . . . , n — 1, 1 < j < n, and r = J2 a a o-o~- 
Then 

[t with n inserted into fs part] = a a [a with n inserted into j's part]. 



One more lemma is quite useful for computations. 
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Lemma 2.4 ( [KW061 Lemma 4.1]). If a planar partition o contains only singleton and 
doubleton parts, and a' is the partition obtained from o by deleting all the singleton parts, 
then the rows of the matrices V^> for a and a' are equal, in the sense that they have the same 
non-zero entries (when the columns are matched accordingly by deleting the corresponding 
singletons) . 

The above lemmas can be used to recursively rewrite a non-planar partition r as an 
equivalent linear combination of planar partitions. As a simple example, to reduce 13 1 245, 
start with the equation from ILemma 2.11 and, using ILemma 2.3[ adjoin a 5 to every part 
containing 4, yielding 

13|245 = 1|2345 + 2|1345 + 3|1245 + 45|123 - 12|345 - 145|23. 

2.4 Connection probabilities 

For a partition r on 1, . . . , n we define 

l = e n a) 

F {i,j}£F 

where the sum is over those spanning forests F of the complete graph on n vertices 1, . . . , n 
for which trees of F span the parts of r. 

This definition makes sense whether or not the partition r is planar. For example, 

-^1|234 = -^2,3-^3,4.+ -^2,3-^2,4 + -^2,4-^3,4 and L 13 | 2 4 = Li^L 2 ,4. 

Recall that Pr(cr) = Pr(cr)/ Pr (uncrossing). 
Theorem 2.5 (Theorem 1.2 of [KW06j ). 

Pr(a)=^pWL T . 

T 

3 Tripartite pairing partitions 

Recall that a tripartite partition is defined by three circularly contiguous sets of nodes R, 
G, and B, which represent the red nodes, green nodes, and blue nodes (a node may be 
split between two color classes), and the number of nodes of the different colors satisfy the 
triangle inequality. In this section we deal with tripartite partitions in which all the parts 
are either doubletons or singletons. (We deal with tripod partitions in the next section.) By 
ILemma 2.41 above, in fact additional singleton nodes could be inserted into the partition at 
arbitrary locations, and the L-polynomial for the partition would remain unchanged. Thus 
we lose no generality in assuming that there are no singleton parts in the partition, so that 
it is a tripartite pairing partition. This assumption is equivalent to assuming that each 
node has only one color. 
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Theorem 3.1. Let a be the tripartite pairing partition defined by circularly contiguous sets 
of nodes R, G, and B, where \R\, \G\, and \B\ satisfy the triangle inequality. Then 








Lr,g 


Lr,b 


Pr(cx) = Pf 


—Lg,r 





Lg,b 




—Lb,r 








Here Lb,g is the submatrix of L whose columns are the red nodes and rows are the green 
nodes. Similarly for Lrb and Lq,b- Also recall that the Pfaffian Pf M of an antisymmetric 
In x In matrix M is a square root of the determinant of M, and is a polynomial in the 
matrix entries: 



Pf M 



) -^^7ri,7r2-^^7r3,7r4 



• • M. 



T2n-l,W2n 



±Vdet M, 



(2) 



permutations ir 

7T1 <7T2 , ■ ■ ■ > 7I "2n- 1 <71"2n 

7ri<7r3<---<7r2„_i 



where the sum can be interpreted as a sum over pairings of {1, . . . ,2n}, since any of the 
2 n n\ permutations associated with a pairing {{n~i, 7r 2 }, . . . , {7r 2n _i, ^2n}} would give the same 
summand. 



In Appendix B there is a corresponding formula for tripartite pairings in terms of the 
matrix R of pairwise resistances between the nodes. 

Observe that we may renumber the nodes while preserving their cyclic order, and the 
above Pfaffian remains unchanged: if we move the last row and column to the front, the sign 
of the Pfaffian changes, and then if we negate the (new) first row and column so that the 
entries above the diagonal are non-negative, the Pfaffian changes sign again. 

As an illustration of the theorem, we have 




Pr(16|23|45) =Pf 











Ll,3 


^1,4 


Ll t 5 


Lifi 










-^2,3 


^2,4 


-^2,5 


L-2,6 


Pf 


— -^1,3 


— -^2,3 








£3,5 


-^3,6 


— -^1,4 


— -^2,4 








-^4,5 






— -^1,5 


— -^2,5 _ 


'£3,5 


— -^4,5 










_ — -^1,6 


— L 2 fi - 


'£3,6 


— £4,6 








-^1,3 


^2,5-^4,6 




5-^4,5 


— L 14 L 2i 5 


-^3,6 + L 


1,4-^: 



(3) 



-^1,5-^2,3-^4,6 + -^1,5-^2,4-^3,6 + -^1,6-^2,3-^4,! 



'3,5 

-^1,6-^2,4-^3,5- 



Note that when one of the colors (say blue) is absent, the Pfaffian becomes a determinant 
(in which the order of the green vertices is reversed). This bipartite determinant special case 
was proved by Curtis, Ingerman, and Morrow [CIM981 Lemma 4.1] and Fomin |Fom01t 
Eqn. 4.4]. See 



8 for a (different) generalization of this determinant special case. 



Proof of \Theorem 3.1\ From Theorem 2.51 we are interested in computing the non-planar 
partitions r (columns of V®) for which Vajr 7^ 0. 
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When we project r, if r has singleton parts, its image must consist of planar partitions 
having those same singleton parts, by the lemmas above: all the transformations preserve 
the singleton parts. Since o consists of only doubleton parts, because of the on the number 
of parts, Va)r is non-zero only when r contains only doubleton parts. Thus in ILemma 2.11 
we may use the abbreviated transformation rule 

13|24 -> -14|23 - 12|34. (4) 

Notice that if we take any crossing pair of indices, and apply this rule to it, each of the 
two resulting partitions has fewer crossing pairs than the original partition, so repeated 
application of this rule is sufficient to express r as a linear combination of planar partitions. 

If a non-planar partition r contains a monochromatic part, and we apply Rule (jlj) to it, 
then because the colors are contiguous, three of the above vertices are of the same color, 
so both of the resulting partitions contain a monochromatic part. When doing the trans- 
formations, once there is a monochromatic doubleton, there always will be one, and since a 
contains no such monochromatic doubletons, we may restrict attention to columns r with 
no monochromatic doubletons. 

When applying Rule (j3]) since there are only three colors, some color must appear twice. 
In one of the resulting partitions there must be a monochromatic doubleton, and we may 
disregard this partition since it will contribute 0. This allows us to further abbreviate the 
uncrossing transformation rule: 

redi x \ red2 y — > — redi y \ red2 x, 

and similarly for green and blue. Thus for any partition r with only doubleton parts, none 
of which are monochromatic, we have Va)r = ±1, and otherwise Vajr = 0. 
If we consider the Pfaffian of the matrix 






Lr,g 


Lr,b 


—Lg,r 





Lg,b 


—Lb,r 


—Lb,g 






each monomial corresponds to a monomial in the L-polynomial of a, up to a possible sign 
change that may depend on the term. 

Suppose that the nodes are numbered from 1 to 2n starting with the red ones, continuing 
with the green ones, and finishing with the blue ones. Let us draw the linear chord diagram 
corresponding to a. Pick any chord, and move one of its endpoints to be adjacent to its 
partner, while maintaining their relative order. Because the chord diagram is non-crossing, 
when doing the move, an integer number of chords are traversed, so an even number of 
transpositions are performed. We can continue this process until the items in each part of 
the partition are adjacent and in sorted order, and the resulting permutation will have even 
sign. Thus in the above Pfaffian, the term corresponding to a has positive sign, i.e., the 
same sign as the a monomial in cr's L-polynomial. 
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Next we consider other pairings r, and show by induction on the number of transpositions 
required to transform r into cr, that the sign of the r monomial in cr's L-polynomial equals 
the sign of the r monomial in the Pfaffian. Suppose that we do a swap on r to obtain a 
pairing r' closer to a. In cr's L polynomial, r and r' have opposite sign. Next we compare 
their signs in the Pfaffian. In the parts in which the swap was performed, there is at least one 
duplicated color (possibly two duplicated colors). If we implement the swap by transposing 
the items of the same color, then the items in each part remain in sorted order, and the sign 
of the permutation has changed, so r and r' have opposite signs in the Pfaffian. 

Thus cr's L-polynomial is the Pfaffian of the above matrix. □ 



4 Tripod partitions 



In this section we show how to compute Pr (cr) for tripod partitions cr, i.e., tripartite parti- 



tions a in which one of the parts has size three. The three lower-left panels of Figure 2 and 



the left panels of Figure 6 and Figure 7 show some examples. 



4.1 Via dual graph and inverse response matrix 

For every tripod partition cr, the dual partition cr* is also tripartite, and contains no part of 
size three. As a consequence, we can compute the probability Pr (cr) when a is a tripod in 
terms of a Pfaffian in the entries of the response matrix L* of the dual graph Q*: 

Pr(crin^) Pr(cr*in£*) Pr(12 - • • n in Q) 



Pr cr 



Pr(l|2| • • • \n in G) Pr(l|2| • • • \n in Q*) Pr(l|2| • • • \n in Q) ' 



The last ratio in the right is known to be an (n — 1) x (n — 1) minor of L (see e.g., § 8); it 
remains to express the matrix L* in terms of L. 

Let i' be the node of the dual graph which is located between the nodes i and i + 1 of Q . 

Lemma 4.1. The entries of L* are related to the entries of L as follows: 

L*',f = - ^+i)L _1 (5j - <5j+i). 

Here even though L is not invertible, the vector 5j — 5j + \ is in the image of L and 5i — 5i+i 
is perpendicular to the kernel of L, so the above expression is well defined. 

Proof. From |K\\'(l(i. Proposition 2.9], we have 

where Rij is the resistance between nodes i and j. From [KWQ61 Proposition A. 2], 

The result follows. □ 
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4.2 Via Pfaffianiod 



In § 4.1 we saw how to compute Pr(a) for tripartite partitions er. It is clear that the formula 
given there is a rational function of the L^s, but from lTheorem 2."5l we know that it is in 
fact a polynomial in the L^'s. Here we give the polynomial. 

We saw in §_3 that the Pfaffian was relevant to tripartite pairing partitions, and that this 
was in part because the Pfaffian is expressible as a sum over pairings. For tripod partitions 
(without singleton parts), the relevant matrix operator resembles a Pfaffian, except that it is 
expressible as a sum over near-pairings, where one of the parts has size 3, and the remaining 
parts have size 2. We call this operator the PfafRanoid, and abbreviate it Pfd. Analogous 
to (j2J), the Pfaffianoid of an antisymmetric (2n + 1) x (2n + 1) matrix M is defined by 



PfdM 



E 



-2 >< M *2n 



(5) 



permutations n 

TL < 7r 2 , ■ ■ ■ ,T2n - 3 < 7I "2n - 2 
1"2n-l<7r2n+l 
7ri<7r3<---<7T2n-3 



where the sum can (almost) be interpreted as a sum over near-pairings (one tripleton and rest 
doubletons) of {1, . . . , In + 1}, since for any permutation associated with the near-pairing 
{{tti, 7r 2 }, . . . , {vr2 n _3, TC 2 n-2} , {^n-i, ^n, 7i~2n+i}}, the summand only depends on the order 
of the items in the tripleton part. 

The sum-over-pairings formula for the Pfaffian is fine as a definition, but there are more 
computationally efficient ways (such as Gaussian elimination) to compute the Pfaffian. Like- 
wise, there are more efficient ways to compute the Pfaffianoid than the above sum-over-near- 
pairings formula. For example, we can write 



PfdM 



E 



\a+b+c 



(M a , b M biC + M h , c M a ,c + M a:C M a , b ) Pf [M s {a, b, c}}, (6) 



l<a<b<c<2n+l 



where M \ {a, b, c} denotes the matrix M with rows and columns a, b, and c deleted. It is 
also possible to represent the Pfaffianoid as a double-sum of Pfaffians. 

The tripod probabilities can written as a Pfaffianoid in the L^-'s as follows: 

Theorem 4.2. Let a be the tripod partition without singletons defined by circularly con- 
tiguous sets of nodes R, G, and B, where \R\, \G\, and \B\ satisfy the triangle inequality. 
Then 



Pi (a) 



\sum of items in it's tripleton part 



Pfd 






Lr,g 


Lr,b 


—Lg,r 





Lg,b 


—Lb,r 


—Lb,g 






The proof of lTheorem 4.21 is similar in nature to the proof of lTheorem 3.11 but there are 



more cases, so we give the proof in |Appendix A 



Unlike the situation for tripartite partitions, here we cannot appeal to ILemma 2.41 to 
eliminate singleton parts from a tripod partition, since ILemma 2.41 does not apply when 
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there is a part with more than two nodes. However, any nodes in singleton parts of the 
partition can be split into two monochromatic nodes of different color, one of which is a leaf. 
The response matrix of the enlarged graph is essentially the same as the response matrix of 
the original graph, with some extra rows and columns for the leaves which are mostly zeroes. 
ITheorem 4.21 may then be applied to this enlarged graph to compute Pr (<r) for the original 
graph. 



5 Irreducibility 

Theorem 5.1. For any non- crossing partition a, Pr(a) is an irreducible polynomial in the 
L{j s. 

By looking at the dual graph, it is a straightforward consequence of ITheorem 5.11 that 
Pr(er)/Pr(12 • • -n) is an irreducible polynomial on the pairwise resistances. In contrast, 
for the double-dimer model, the polynomials Pr(a) sometimes factor (the first, second, and 



fourth examples in § 6 factor). 



Proof of \Theorem 5.11 Suppose that Pr (a) factors into Pr (a) = P\P 2 where Pi and P2 are 
polynomials in the -Z^j's. Because Pr(cx) = ^ t P£tPt and each L T is multilinear in the 
Ljj's, we see that no variable occurs in both polynomials Pi and Pi- 

Suppose that for distinct vertices i,j,k, the variables L^j and both occur in Pr(cx), 
but occur in different factors, say Lij occurs in Pi while occurs in P 2 . Then the product 
Pr(cr) contains monomials divisible by L^jL^f.- If we consider one such monomial, then the 
connected components (with edges given by the indices of the variables of the monomial) 
define a partition r for which V&jr 7^ and for which r contains a part containing at least 
three distinct items i, j, and k. Then L T contains Lj^, so Lj^ also occurs in one of Pi or 
P2, say (w.l.o.g.) that it occurs in Pi. Because L T contains monomials divisible by L^jLj^, 
so does Pr(er), and hence Pi must contain monomials divisible by L$jZy&. But then P1P2 
would contain monomials divisible by LijLi^Lj^, but Pr(cr) contains no such monomials, a 
contradiction, so in fact Lij and L^f. must occur in the same factor of Pr (er). 

If we consider the graph which has an edge {i, j} for each variable Li j of Pr(cr), we aim 
to show that the graph is connected except possibly for isolated vertices; it will then follow 
that Pr(cr) is irreducible. 

We say that two parts Q\ and Q 2 of a non-crossing partition a are mergeable if the 
partition a \ {Qi, Q 2 } U {Qi UQ 2 } is non-crossing. It suffices, to complete the proof, to show 
that if Qi and Q 2 are mergeable parts of a, then Pr(cr) contains L a c for some a G Q\ and 
c G Q 2 . 

Suppose Qi and Q 2 are mergeable parts of a that both have at least two items. When the 
items are listed in cyclic order, say that a is the last item of Q\ before Q 2 , bis the first item of 
Q 2 after Qi, c is the last item of Q 2 before Qi, and d is the first item of Q\ after Q 2 . Let r be 
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the partition formed from a by swapping c and d. Let a* = a\{Qi, Q2}, and let A = Qi\{d} 
and B = Q 2 \ {c}. Then a = a* U {A U {d}, B U {c}} and r = a* U {AU {c}, B U {d}}. Then 

r — > — o" — (cr* U {A U 5, {c, 4}) 

+ ((T*U{AU5U {rf}, {c}}) + ((7*U{ J 4UBU {c}, {4}) 

+ ((7* U{A,BU {c, 4}) + (a* U {A U {c, 4, 5}). 

Each of the partitions on the right-hand side is non-crossing, so Va)r = —1, so in particular 
L a c occurs in Pr(er). 

Now suppose that a contains a singleton part {a} and another part Q2 containing at 
least three items b, c, d, where b, c, and d are the first, second, and last items of the part Q2 
as viewed from item a. Let a* = a \ {{a}, Q2} and C = Q 2 \ {b, d}. Let r be the partition 

T = a*U{{a}UC,{b,d}}. 

Now 



r -> a + (a* U {{b}, {a, d} U C}) + (a* U {{4, {a, b} U C}) + (a* U {C, {a, 6, 4}) 

- (a* U {{&} U C, {a, 4}) - {a* U {{a, b}, {d} U C}). 

The second, third, fourth, fifth, and sixth terms on the RHS contribute nothing to 
because their restrictions to the intervals [b,b], [d,d], (b,d), [b,d), and (6, d] respectively are 
planar and do not agree with a. Thus Va)r = 1, and hence L a c occurs in Pr(cr). 

Finally, if a contains singleton parts but no parts with at least three items, then Pr(cr) 
is formally identical to the polynomial Pr(a*) where a* is the partition with the singleton 
parts removed from a, and we have already shown above that the polynomial Pr(cx*) is 
irreducible. □ 



6 Tripartite pairings in the double- dimer model 

In this section we prove a determinant formula for the tripartite pairing in the double-dimer 
model. 

Theorem 6.1. Suppose that the nodes are contiguously colored red, green, and blue (a color 
may occur zero times ), and that a is the ( unique ) planar pairing in which like colors are not 
paired together. Let o~i denote the item that a pairs with item i. We have 

Pr(o") = detfl^ j colored diSerently^ijjj^^as,...,^-!' 
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For example, 




Pr(ilS) 



x 1A 

X 3 , 2 



(this first example formula is essentially Theorems 2.1 and 2.3 of [Kuo04] . see also |Kuo06] 
for a generalization different from the one considered here) 




Pn 



1 I 3 I 5 

2 I 6 I 4 



-^1,2 





^1,4 





-^3,6 





-^5,2 





-^5,4 




Pn 



1 I 3 I 5^ 

2 I 4 I 6; 



Pr( 



1 


3 


5 1 


2 


8 


6 1 



-i^l,2 


^1,4 










-^3,4 


-^3,6 




^5,2 





-^5,6 




-^1,2 








^1,4 





-^3,8 


-^3,6 








-^5,8 


-^5,6 





^7,2 








-^7,4 


X\,2 





-^1,4 


-^1,6 





^3,8 





-^3,6 


-^5,2 


^5,8 


-^5,4 





^7,2 





-^7,4 


-^7,6 



PrT 1 1 3 I 5 I 7 ' 

rl l 2 I 8 I 4 I 6 , 



Proof. Recall our theorem from [KW06j . which shows how to compute the "X" polynomials 
for the double-dimer model in terms of the "L" polynomials for groves: 

Theorem 6.2 (Kenyon- Wilson '06). If a partition a contains only doubleton parts, then if 
we make the following substitutions to the grove partition polynomial Pr(cr): 







if i and j have the same parity 
otherwise 



(_l)(l*-i|-l)/2 Xi 

then the result is the double-dimer pairing polynomial Pr (a) , when we interpret a as a pairing. 
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Thus our Pfaffian formula for tripartite groves in terms of the £j,/s immediately gives a 
Pfaffian formula for the double-dimer model. For the double-dimer tripartite formula there 
are node parities as well as colors (recall that the graph is bipartite). Rather than take a 
Pfaffian of the full matrix, we can take the determinant of the odd/even submatrix, whose 
rows are indexed by red-even, green-even, and blue-even vertices, and whose columns are 
indexed by red-odd, green-odd, and blue-odd vertices. For example, when computing the 
probability 

rl I 8 I 4 I 6 I 2 )i 

nodes 1, 2, and 3 are red, 4 and 5 are green, and 6, 7, and 8 are blue; the odd nodes are 
black, and the even ones are white. The L-polynomial is 












£l,4 


£1,5 


£l,6 


£1,7 


£1,8 











£2,4 


£2,5 


£2,6 


£2,7 


£2,8 











£3,4 


£3,5 


£3,6 


£3,7 


£3,8 


— £4,1 


-£4,2 


-£4,3 








£4,6 


£4,7 


£4,8 


— £5,1 


— £5,2 


~ £5,3 








£5,6 


£5,7 


£5,8 


— -^6,1 


— -^6,2 


— Lq,3 


~ -^6,4 


— ^6,5 











— £7,1 


"-^7,2 


— £7,3 


— £7,4 


— £7,5 











— £8,1 


—Ls,2 


— -^8,3 




— £8,5 












Next we do the substitution Ljj — > when i + j is even, and reorder the rows and columns 
so that the odd nodes are listed first. Each time we swap a pair of rows and do the same 
swap on the corresponding pair of columns, the sign of the Pfaffian changes by —1. Since 
there are 2n nodes the number of swaps is n{n — l)/2. If n is congruent to or 1 mod 4 the 
sign does not change, and otherwise it does change. After the rows and columns have been 
sorted by their parity, the matrix has the form 

±L ,e 
TL e ,o J ' 

where O represents the odd nodes and E the even nodes, and where the individual signs 
are + if the odd node has smaller index than the even node, and — otherwise. The Pfaffian 
of this matrix is just the determinant of the upper-right submatrix, times the sign of the 
permutation 1, n + 1, 2, n + 2, . . . , n, 2n, which is (-l)™^- 1 )/ 2 . This sign cancels the above 
(— 1)"(" -1 )/ 2 sign. In this example we get 



det 






£1,4 


£1,6 


£1, 





£3,4 


£3,6 


£3, 


-£5,2 





£5,6 


£5, 


—£7,2 


—£7,4 









Next we do the L i: j — ► (— l)d* ^ ^^JQj substitution. The i,j entry of this matrix is 
(— iy >j Lij. Each time that % or j are incremented or decremented by 2, the (— l)(l l ~^V 1 )/ 2 
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sign will flip, unless the (— 1) %> ^ sign also flips. After the substitution, the signs of the Xij 
are staggered in a checkerboard pattern. If we then multiply every other row by —1 and 
every other column by —1, the determinant is unchanged and all the signs are +. In the 
example we get 



det 






—Xi^ 


-^1,6 









X\£ 


-^1,6 







-^3,4 


~^3,6 


-^3,8 


= det 





-^3,4 


-^3,6 


X3, 







-^5,6 




-^5,2 





^5,6 




^7,2 


^7,4 










_^7,2 


^7,4 









There is then a global sign of (— l) a where the sign of the pairing a is the sign of sign of the 
permutation of the even elements when the parts are arranged in increasing order of their 
odd parts. In our example, the sign of g 1 4 1 1 1 2 * s the s ^S n °f 8462, which is —1. This global 
sign may be canceled by reordering the columns in this order, i.e., so that the pairing a can 
be read in the indices along the diagonal of the matrix, which for our example is 



det 



X\ t s 


Xi } 4 


X\fi 





-^3,8 


^3,4 


-^3,6 





^5,8 





^5,6 


-^5,2 





^7,4 





Xr,2 



□ 



7 Reconstruction on the "standard graphs" S ra 

Given a planar resistor network, can we determine (or "reconstruct") the conductances on 
the edges from boundary measurements, that is, from the entries in the L matrix? 

While reconstruction is not possible in general, each planar graph is equivalent, through 
a sequence of electrical transformations, to a graph on which generically the conductances 



can be reconstructed. Let S n denote the standard graph on n nodes, illustrated in Figure 5 
for n up to 6. Every circular planar graph with n nodes is electrically equivalent to a minor 
of a standard graph E n . 

Here we will use the Pfaffian formulas to give explicit formulas for reconstruction on 
standard graphs. For minors of standard graphs, the conductances can be computed by 
taking limits of the formulas for standard graphs. 

Curtis, Ingerman and Morrow |CIM98j gave a recursive construction to compute con- 
ductances for standard graphs from the L-matrix. Card and Muranaka |CM02j give another 
way. Russell [Rus03j shows how to recover the conductances, and shows that they are ratio- 
nal functions of L-matrix entries. However the solution is sometimes given parametrically, 
as a solution to polynomial constraints, even when graph is recoverable. 

For a vertex v G S n we define tt v to be the tripod partition of the nodes indicated in 



Figure 6, with a single triple connection joining the nodes ?;_> horizontally across from v and 
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the two nodes u-f , vertically located from v (in the same column as v), and the remaining 
nodes joined in nested pairs between and v^, and v\, and v\ and (with up to two 
singletons if t>^, and/or t>| have an odd number of nodes between them). 

Similarly, for a bounded face / of E n define tij to be the tripartite partition of the nodes 



indicated in Figure 6 , It has three nested sequences of pairwise connections (with two of the 
nested sequences going to the NE and SE, possibly terminating in singletons). We think of 
the unbounded face as containing many "external faces," each consisting of a unit square 
which is adjacent to an edge of E n . For each of these external faces, we define 7T/ in the same 
manner as for internal faces. For the external faces / on the left of E„, the "left-going" nested 
sequence of 717 is empty. For the other external faces /, the partition 7T/is(l,n|2,n — 1| • • • ), 
independent of /. 



v 1 



/ 



r 



Figure 6: Tripod partition (left) and tripartite partition (right) on the standard graph E n . 

Observe that for the standard graphs E n , there is only one grove of type ir v or of type iif. 
Let a e denote the conductance of edge e in E n . Each Z nv and Z nf is a monomial in these 
conductances a e . To simplify notation we write Z v = Z 7Tv and Zf = Z Wf . 

Each conductance a e can be written in terms of the Z v and Zf. 

Lemma 7.1. For an edge e of the standard graph E n , let v\ and v% be the endpoints of e, 
and let fi and f'2 be the faces bounded by e. We have 

7 7 
a e = — — — — . 

z h z h 

Proof. A straightforward inspection of the various cases. □ 

Combining this lemma with the results of Sections H] and [3X we can write each edge 
conductance rational function in the L^s. Since the Z„'s and Zfs are irreducible 

by lTheorem 5.11 this formula is the simplest rational expression for the a e 's in terms of the 
Ljj's. 
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8 Minors of the response matrix 

We have the following interpretation of the minors of L. 

Theorem 8.1. For general graphs (not necessarily planar), suppose that A, B, C, and D 
are pairwise disjoint sets of nodes such that \A\ = \B\ and AUflUCUD is the set of all 
nodes. Then the determinant of L(auc),{buc) is given by 



det [L it , 



1 i=ai,...,0| A |,c 1 ,...,C| | _ XLl -1 )"" P r [ fl l> ^(1)1 - ' ' \ a \A\, K(\A\)\dl\ ■ ■ ■ \d\ D \ 



:,3ij=b 1 ,...,b m ,c 1 ,...,c ]c \ ~ Pr [uncrossing] 

where the nodes of C may appear in any of the above parts. 



In Appendix B equation ( fT2i) . there is a corresponding formula in terms of the pairwise 
resistances between nodes. 

For an example of the Theorem, if there are 6 nodes, then 

detL (lj2! 3), (3 ,4,5) =Pr(15|24|6) - Pr(14|25|6) 

=Pr(153|24|6) + Pr(15|243|6) + Pr(15|24|63) 
- Pr(143|25|6) - Pr(14|253|6) - Pr(14|25|63), 

which for circular planar graphs is just Pr(15|243|6). 

When C = 0, this determinant formula is equivalent to Lemma 4.1 of Curtis- Ingerman- 
Morrow [CIM98], though their formulation is a bit more complicated. The formula Lij = 
Pr(i, j | rest singletons) [KW06t Proposition 2.8] is a further specialization, with A = {i} and 

" {./)• 



det L(a\jc),(b\jc) 



La,b 


La,c 




Lc,b 


Lc,c 





I 

Lc.b Lac 



Proof of Theorem. We assume first that La,c is non-singular. By standard linear algebra 

La,b — La,cL c x c L cb La,cL c x c 
I 

= det[L AiB - La,cL c ] c Lc,b] det L c , c (7) 

Since this is essentially Schur reduction, La,b — La,cL c x c Lc,b is the A, B submatrix of the 
response matrix when nodes in C are redesignated as internal, so by Lemma 4.1 of Curtis- 
Ingerman- Morrow |CIM 98j. 

x signed sum of pairings from A to B when C is internal 

det [La,b — L AC L C c L c ,b\ = ; : -pr. — — \ • (8) 

uncrossing when C is internal 

If we glue the nodes not in C together, the response matrix of the resulting graph has 
Lc,c as a co-dimension 1 submatrix, so by Lemma A.l of Kenyon- Wilson [KW06J, 

forests rooted at A U B U D uncrossing when C is internal . . 

det L c ,c = ; = ; • (9) 

uncrossing uncrossing 

Combining Equations EJ and [9] gives the result for nonsingular Lq\c- 

The case of singular Lq\c can be obtained as a limit of the above nonsingular case. □ 
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9 Carroll-Speyer groves 

Here we study the groves of Carroll and Speyer. For Carroll and Speyer's work, the relevant 



graph is a triangular portion of the triangular grid, shown in Figure 7 Carroll and Speyer 
computed the number of groves on this graph which form a tripod grove (for N even) or 
a tripartite grove (for N odd). The relevant tripod or tripartite partition is the one for 
which the three sides of the triangular region form the three color classes, and each part 
connects nodes with different colors. For the case N = Q, the relevant tripod partition is 



1, 17|2, 16 1 3, 9, 15|4, 8|5, 7|10, 14| 11, 13 1 6 1 12 1 18 , and an example grove is shown in Figure 7 
The number of such groves turns out to always be a power of 3, specifically, when there are 
3iV nodes, there are 3^ / 4 J groves. In this section we consider these graphs as a test case 
for our methods for counting groves. There is much that we can compute, but we do not 
know how at present to obtain a second derivation of the power-of-3 formula. 

We need the entries of the L matrix in order to compute the connection probabilities using 



the Pfafnan and Pfaffianoid formulas presented in § 3 and § 4 To compute the tripartite 
connection probabilities, we need those entries of the L matrix whose indices come from 
different sides of the triangle. From symmetry considerations, it suffices to consider the 



entries between the first two sides. In the N = 6 example from |Figure 7 the submatrix of 



the L matrix with rows indexed by the nodes on side 1 (excluding corners) and columns 




Figure 7: Carroll-Speyer graph for N = 6 (left) and N = 7 (right), each shown with one of 
its Caroll-Speyer groves. The graphs have n = 3N nodes; for even N the grove partition 
type is a (generalized) tripod, while for odd N the grove partition type is a (generalized) 
tripartite pairing. The number of Carroll-Speyer groves is 3^ / 4 J [CS04]. 
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indexed by the nodes on side 2 (excluding corners) is given by 







-^1,9 


£1,10 


-^1,11 


L<2,7 


-^2,8 


-^2,9 


-^2,10 


-^2,11 


£3,7 


£3,8 


-^3,9 


-^3,10 


-^3,11 


L<4,7 


£4,8 


-^4,9 


-^4,10 


-^4,11 


£5,7 


£5,8 


-^5,9 


£5,10 


-^5,11 



31 


25 


23 


25 


31 


945,6 


2364 
445 


15^6 


2364 
355 


9 i 6 


2364 
87 


9456 
53 


394 
97 


9456 
23 


2364 
23 


1576 
529 


394 
3043 


788 
53 


394 
445 


1576 
25 



2364 9456 
11167 529 

9456 2364 
1 

-31 

361 - 
-2015 
5297 



391 
87 

1576 
-1 

24 
-208 
888 
-2015 



2364 
1 

-16 

81 
-208 
361 



2 lf 4 

9456. 
-1 

8 

-16 
24 
-31 



n -1 



We explain here why the inverse of this submatrix is integer-valued, and how to interpret 
the entries. 

Recall ITheorem 8.11 on minors of the L matrix. Letting A, B, and C denote the nodes 
on the first, second, and third sides respectively, 



Z (nodes of A paired with nodes of B, nodes of C singletons) 



1 



Likewise 



det^,^. 



iGA\{i } 



Z (uncrossing) Z (uncrossing) 

nodes of A \ {io} paired with nodes of B \ {jo}, 
nodes of C U j } singletons 



?es\{io} 2, (uncrossing) 

Thus the %q, jo entry of the inverse of the above matrix is 



Jo,«o 



nodes of A \ {i } paired with nodes of B \ {jo}. 
nodes of C U {i , j } singletons 



When there are edge weights, the io,jo entry of the inverse matrix will be given by the 
corresponding polynomial in the edge weights. 

To get the normalized probability of the tripartite partition (for odd N) , the Pfafnan we 
need is 

L AtB 









La,b 


La,c 


Pf 


—Lb,a 





Lb,c 




—Lc,a 


—Lc,b 






Pf 



-T T 

A,B 



tT ■ 

A,B 



— L 



A,B 



-T T 

^A,B 



La,b 





which in the case N = 7 gives Pr (tripartite grove) = 531441/135418115000. The calculations 
for the tripod partitions for even N is similar, except that we take a Pfaffianoid rather than 
a Pfafnan. 
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To compute the number (as opposed to probability) of groves of a given type, we also 
need the number of spanning forests rooted at the nodes. The number of spanning forests 
may be computed from the graph Laplacian using the matrix-tree theorem, which yields the 
following formula 

[ [ (6-a-a- 1 -/3-/5- 1 -7-7" 1 ) 

{a,P,f} 
a 3JV >=l 
(a//J)*=l 
a/?7=l 
a,/3,7 distinct 

(see [KPWOO, § 6.9] for the derivation of a similar formula). In the case N = 7 this formula 
yields 135418115000, so there are 531441 = 3 12 = 3^ tripartite groves, in agreement with 
Carroll and Speyer's formula. Is it possible to derive the 3^ / 4 J formula using this approach? 



A Pfaffianoid formula for tripod partitions 



Proof of Theorem 4-2 . Any column partition r contributing to a will have (n— 1)/2 parts (as 
a does) and no singleton parts, and as such it will consist of a single tripleton part together 
with doubleton parts. To determine what r contributes to a, we may use the following 
abbreviated rules. For two crossing doubletons, as in (j3J) we use 

13|24^ — 12|34 — 14|23. 

For a crossing doubleton and a tripleton, recall (Lemmas 12.11 and 12. 3j) that 

13|24 = 1|234 + 2|134 + 3|124 + 4|123 - 12|34 - 14|23 
13|245 = 1|2345 + 2|1345 + 3|1245 + 45|123 - 12|345 - 145|23 

so we may use the rule 

13|245 -> 45|123 - 12|345 - 23|154 

After we apply the 13|245 rule, let us consider another doubleton part. If the doubleton 
part did not cross 13 or 245, it will cross none of 45, 123, 12, 345, 23, or 154. Otherwise it 
is one of the following forms: 
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1.5, 3.5 
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1.5, 4.5 
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n 


2.5, 4.5 
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y 


n 


y 
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y 


3.5, 4.5 


n 


y 


y 


n 


n 


y 


n 


y 



In each case the doubleton part crosses at most as many other parts in the new partitions 
as it did in the old partition. Thus the number of crossing parts strictly decreases. 

After we apply the 13 1 24 rule, we saw already that the number of crossings amongst 
doubleton parts strictly decreases. Now let us consider crossings amongst doubleton parts 
and the tripleton part. The tripleton part divides the vertices into three sectors. The 
distribution of 1, 2, 3, 4 amongst these sectors is one of 

• All four endpoints in same sector. Before rule neither of the chords 13,24 cross the 
tripleton, nor do any of them after the rule. 

• Three endpoints in same sector. Before rule one chord crosses tripleton, after rule one 
chord crosses tripleton. 

• At least one sector contains exactly two endpoints. If we apply the rule then it must be 
that the chords crossed, so both chords cross the tripleton. After rule, in one partition 
two chords cross tripleton, in other partition neither chord crosses tripleton. 

In any event, the total number of crossing parts strictly decreases. 

We may repeatedly apply the 13 1 245 rule and 13 1 24 rule, in any order, until no two parts 
cross, and we are guaranteed to obtain a linear combination of planar partitions containing 
a tripleton part and rest doubleton parts. 

Recall the rule 

13|245 -> 45|123 - 12|345 - 23|154. 

Suppose that at least one of these parts has two or more vertices of the same color, say red. 
We have the following possibilities: 

12|345 - 23|154 

r r r r 

12|345 - 23|154 
12|345 - 23|154 
121345 - 231154 



r r r 

13 1 245 

r r 

13 1 245 

r r r 

13|245 

r r r 



r r r 

45|123 

r r 

45|123 

r r r 

451123 



131245 -> 451123 - 
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(the remaining cases contain more reds than these). In each case, each of the resulting 
partitions has a part with two or more vertices of the same color, so by induction, such 
partitions will not contribute to ex. Thus we may restrict attention to partitions r with 
bichromatic doubleton parts and a trichromatic tripleton part. 
We take another look at the rule 

13|245 -> 45|123 - 12|345 - 23|154. 

Two colors occur twice. There are three possibilities (up to renamings of color): 

5,1 red and 3,4 blue in which case we may use the rule 13|245 — > 45| 123 
(*) 5,1 red and 2,3 blue in which case we may use the rule 13 1245 — > — 12 1 345 
(**) 1,2 red and 3,4 blue in which case we may use the rule 13 1245 — > — 23 1 154 

(And similarly for the 13|24 rule, either 13|24 -> -12|34 or 13|24 -> -14|23.) 

Thus each maximally multichromatic partition r contributes either +1 or —1 to o. 

Let us compare the contributions to o of partitions r which contain doubleton parts that 
do not cross one another, and a tripleton part which may cross the doubleton parts. Any 
such partition r may be obtained from a by repeatedly transposing items of the tripleton 
with their neighbors. From the above rules (*) and (**) we see that each such transposition 
changes the sign of r's contribution to a. 

Next we consider what happens if we leave the tripleton part alone and only apply the 
13 1 24 rule. Upon deleting the tripleton part, and recalling our earlier analysis of the tripartite 



pairing ([Theorem 3.1j) . we see that each partition r contributes to a a sign which involves 
moving the tripod to its desired location in a, times the sign in the corresponding Pfaffian 
in which rows and columns of items of tripleton part have been deleted. This gives us the 
tripod Pfaffianoid formula: 

/ |jsum of items in cr's tripleton part ^ 

(-l) a+b+c (L atb L btC + L a , b L a , c + L a , c L biC ) Pf(i? s {a}; G \ {b}; B \ {c}). 

aeR,beG,ceB 

□ 



B Tripartite pairings in terms of Pfaffians in the resis- 
tances 

There is a formula analogous to the one in ITheorem 3.11 for the probability of a tripartite 
partition, expressing it in terms of the pairwise electrical resistances Ri j rather than the 
response matrix entries L it j. The normalization is also different, rather than normalizing by 
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Pr(l|2| ■ • • |n), we normalize by Pr(12 
following example: 

Pr(16|23|45) 



•n 



The corresponding formula is illustrated by the 



Pr(123456) 



[f]Pf 






— t + |-Rl,3 
—t + 

— < + 5 #1,5 
— t + 9-^1,6 







-t + 2^2,3 
-t + 2-^2,4 
-t + 2-^2,5 
-t + 0-^2,6 



t — 9-^1,3 

2 : 3 



t - §# 






-t + ~# 
-t + \R. 



3,5 
3,6 



i — 2" -Ri,4 

t — 2-^2,4 




-£ + ^R. 
-t + 2-R. 



4,5 
4,6 



t 


- 2^1,5 


t 


— 2^1>6 


t 


- 2-^2,5 


t 


~~ 2^ 2 > 6 


t 


- 2-^3,5 


t 


~~ 2^3,6 


t 


— 2^ 4 > 5 


t 


— 2"^ 4 > 6 























(10) 



Here the Pfaffian is a polynomial in t, in fact a linear function of t, and the coefficient of the 
linear term gives Pr(16|23|45)/ Pr(123456). 

We prove here this formula when one of the color classes is empty, so that the Pfaffian 
is actually a determinant (see lTheorem B.3l . The general tripartite Pfaffian formula follows 
from the bipartite determinant special case, together with a result that we prove in our 
next article that for any planar pairing, the L-polynomial is a linear combination of such 
determinants [KW08J. 

Recall that any codimension-1 minor of L is the ratio of the spanning trees to the un- 
crossing. Let L be the (n — 1) x {n — 1) submatrix of L obtained by deleting the last row 
and column. Recall that i?y = L^ 1 + Ljj — 2L^j, where L~j is interpreted to be if either 
i = n or j = n (see [KW06] Section A.2]). 

Lemma B.l. For each row of I + LR/2, the row's entries are all the same. 
Proof. Suppose i ^ n. Then 

3 3 

= Li t j[Rj >n + 0] — 25i ; k- 



□ 



So for any % ^ n, the i row of / + LR/2 is constant. By re- indexing the rows, we see that 
this must hold for the n th row as well. 

Lemma B.2. The diagonal entries of RLR are all the same. 
Proof. We have 

y, Ri,jL j:k Rk,i = /jLj} + Ljj - 2L^j)L jik (Li t l + L^\ - 2L~l). 

3,h j,k 
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The third factor contains no j subscripts, and neither does the L~\ term of the first factor, 
so summing over j gives 0, and we may drop the L^l term in the first factor. Similarly, we 
may drop the L~\ term in the third factor. 

2j RijL jtk R k>i = 22(Ljj - 2Lrj)L j<k (L k j c - 2Lj£) 

j,k j,k 

= LjjLj^L^ - 2 ^2 ^j,j L j,kLi t l - 2} u L~jL jik (L^ k - 2Lr^) 

j,k j,k j,k 

= S Lj,j L j,kL k l ~ 2 ^2 L J-j' S '-' ~ 2 ^2 W^k,* ~ 2 ^i,k) 
j,k j k 

= ^l] L 3^k]k 
j,k 

which in particular does not depend upon i: the diagonal terms of RLR are all equal. □ 

Suppose that we add an extra (n + l) st node to the graph where the conductance between 
nodes j and n + 1 is e{l + \ Y^ fc Lj k Rj jk ), where e ~ 0. (Some of these conductances may 
be negative.) From [Lemma B.TI we have 



$ij + 9 L lkRk,i - 1 + - ^2 L i,kRj,k 



h,3 + 9 L 3,kRk,i - ^ 1 + 9 L 3,kRj,k 



so adding up the conductances of these new edges gives e. If we set node n + 1 to be at 1 
volt and node i to be at volts, then to first order, each of the nodes 1, . . . , n be nearly at 
volts. Then the current flowing from node n + 1 to node j is (1 + o(l))e(8 ii j + ~ V^ fc Lj tk Rk,i), 
and the total current from node n + 1 is (1 + o(l))e. Let v denote the voltages of the first 
n nodes. Now we view the first n nodes as a circuit to which we apply the voltages v. The 
current flowing into node j 7^ i is (1 + o(l))| J2 k Lj jk R k ,i, and the current flowing into node i 
is (1 + o(l))e(l + \ J2 k L ijk R kji ) — (1 + o(l))e (the first term is the current from the (n + l) st 
node to node i, and the second term is the current flowing out of node i, i.e., the total current 
flowing in from node (n + 1)). To summarize, we have 



Lv = (l + o(l))-eLR6i, 

where 6i is the i th basis vector. We would like to cancel the L's from this equation, but the 
equation only determines v up to an additive constant: v = \eR5i + o(e) + const. But t>j = 
and R iti = 0, so the additive constant is zero. Knowing the voltages Vj = \sRij +o(e) allows 
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us to compute the currents to second order: the total current from node n + 1 is 



+ 2 ^ L hk R k,i (l - ~Ri,j + o(e)^ - e - — ^2 Ri,j L j,kRk,i + o(e 2 ) 



e 2 

= e-—Q + o(e 2 ) 

where Q is the diagonal entry of RLR. Let R denote the resistances in the extended graph; 
we have 

Ri,n+i = l/e + Q/4 + o(l). 

Let L be the L-matrix of the enlarged graph, L' be the submatrix obtained by deleting 
the (n + l) st vertex, and L'^ 1 be the inverse of L' (sometimes extended to have an all-0 
(n + l) st row and column). The ratio of spanning trees to the uncrossing is 

1 11 

lim - det(— L') = lim 

e-o e £ -oedet(-L'- 1 ) 



now recall {-V ^ = {Ri, n +\ + Rj,n+i ~ Ri,j)/2, so 



= lim -- 



1 



-o e det[l/e + Q/4 - R hJ /2 + o(l)]jzi',V.^' 
This determinant is a linear function of t = 1/e + Q/4, so we may rewrite this limit as 

1 



[t] det[t - ii^-Eli' 



i=l,. ..,n 



Let A and 5 be two subsets of the nodes. We have 



detMti=Kmdet[Ly5li 



limdetL'detfL^ 1 ]^ x (-l)E ae A«+E bes ' 



= lim 



__ det[-l/e - Q/4 + R^/2 + o(l)}^ 



X (_l)EaeA a +Ef> e B b 

[t]det[-t + |^g c 
>]det[i-|i^]^ 



(11) 



where in the above, the rows and columns in A, B, A c , and B c are arranged in sorted order. 
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Note that (TTTT) allows us to rewrite any minor of the L matrix in terms of the pairwise 
resistances between the nodes. Since the denominator of the right-hand side of ( TTTT) is 
Z (uncrossing) jZ (tree), we have 



Z (uncrossing) 



det[Ly]}H = [t] det[-t + \RJS X (-1)£^+E> 6B >. (12 ) 



Z(tree) 

If A = B c , then (fl2|) simplifies to yield 

Theorem B.3. //A and B are disjoint and equinumerous sets of nodes and AU B is the 

set of all nodes, then 

det[L,4H = [t] drt[t - iiy'Ii 
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